Rectangular microwave resonators with magnetic anisotropy. 
Mapping onto pseudo integrable rhombus 



Evgeny N.Bulgakov ^ and Almas F. Sadreev^'^ * 
1 ) Kirensky Institute of Physics, 660036, Krasnoyarsk, Russia 
2) Department of Physics and Measurement Technology, Linkoping University, S-581 83 

Linkoping, Sweden 

o 
o 

(N 

Abstract 

The rectangular microwave resonator filed by a ferrite with uniaxial magnetic 
anisotropy is considered. It is shown that this task can be reduced to an 
empty rhombus resonator with vertex angle defined by external magnetic field 
provided the magnetic anisotropy of the ferrite is strong. Therefore statistics 
of eigen frequencies for TM modes is described by the Brody or semi-Poisson 
distribution with some exceptional cases. 
One of the main research lines in quantum chaos is to investigate the statistics of energy 
levels and eigen functions of quantum systems whose classical counterpart is chaotic. A very 
popular class of systems are two-dimensional Euclidean billiards, which is described by the 
Helmholtz equation 

O " 

O; -V'^{x,y) = Xij{x,y) (1) 

O ■ 

with Dirichlet boundary conditions ipi^y Z/) = fo^' u) ^'^ the boundary of billiard. It was 
O I shown that the eigen value statistics of A obey the statistics of random matrix ensemble 
I 001 • '^^^ distribution function for the amplitudes of the eigenfunctions ip is perfectly well 

described by a Gaussian distribution P,^. Correspondingly the square p = IV'P is describing 
by the well known Porter-Thomas (P-T) distribution |Q 
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P(p) = -=exp(-p/2). (2) 



a • V27rp 

For integrable billiards, for example, rectangular one, the eigen value statistics is de- 
scribing by the Poisson distribution. 

The distribution function of p is given by formula Q 

g{p) = ^K{l-p), 0<p<l (3) 
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where K{m) is the elhptic integral of the first kind and the eigen state is il){x,y) = 
sm{kxx) sin{kyy). The distribution (§) has a stepwise behavior at p = 1 with a value of 
step equaled i^(0)/7r^ resulted by that p{x,y) has hte same maxima. 

In the two decades a subject of these distributions was dominated by theory and com- 
puter simulations. About ten years ago experimentalists found effective techniques to study 
billiards with the help similarly shaped electromagnetic resonators There are two 

ways to obtain precisely the Helmholtz equation (|ip [Q. For the TM modes we have 
ip{x,y) = Ez{x,y), = and A = (w/c)^ — {Trn/d)'^,n = 0,1,2,... where the z-th compo- 
nent of electric field obeys the Dirichlet boundary condition, uj is the angular frequency 
and c is the light velocity, ci is a thickness of the resonator. For the TE modes we have 
y) = Bz{x, y),Ez = but with the Neumann boundary conditions for ip{x, y). In what 
follows we consider here only the TM modes. 

We consider more general case of the Helmgoltz equation 

{/ixx^ + ^^yy^ + (/ixj; + f^y-)-Q^ + ^l^l^^,!/) = 0. (4) 

with the Dirichlet boundary conditions where components of tensor fia,i3 will be defined 
below. Because of term with mixed partial derivatives in @ even for the rectangular billiard 
the solution of this equation can not be presented as ifj^x^y) = ilj{x)'^{y) and, therefore, 
even the rectangular billiard is becoming non integrable one. Below we will show that 
Eq.(^ for the rectangular billiard can be transformed to the isotropic Helmholtz equation 
(H) but with the Dirichlet boundary conditions at the rhombic billiard. Eq.(P describes 
microwave resonator filed by anisotropic magnet. Using of magnets in microwave resonators 
was described already as a way to violate time reversal symmetry Here we consider 



a complete filling of resonator to have homogeneous case for which time-reversal symmetry 
takes place 

Let us write the Maxwell equations for TE modes of electromagnetic fields in two- 
dimensional resonators 

V ■ B = 0, 

V X nEz = -ikB, 
V X H = iknEz, 

B = pH, (5) 

where n is unit vector parallel to electric field directed , H is the magnetic field, B is the 
magnetic induction, k = uj/c and is a eigen frequency with wave number k. In what 
follows we imply the following magnetic anisotropic permeability 

/i = 1 + X, 

/ ~X.xx Xxy 
Xyx Xyy I . (6) 



X 



V xz 



Following to 0] we consider the TM modes of the Maxwell equations (j^). Substituting 
the permeability into the last Maxwell equation we have 
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-^a; \ ^ J_ / fJ-yy fJ-xy \ f k dy 

HyJ D \liyx f^XxJ \ 



(7) 



where 

D ^^xx^^yy f^xy^yx- 

From this equation one can find from the Maxwell equations (|^) the wave equation for 
electric field with eigen values equaled to A = Dk^. 
By the coordinate transformation 



y / \ i-^xy~\~ ^yx 

'^fJ'XX 




(8) 



we can eliminate the cross derivatives in Eq. (^) and reduce the Helmholtz equation to the 
following on 



V'E, + fik'E, = (9) 



where 



^ ^ ~( Z vJa' ^ 

fj'xx^^yy [f^xy + f^yx) /4 

The transformation (H) transforms the rectangle into a particular case of polygon as a 
parallelogram as shown in Fig. |1| (a). 

Below we consider ferrite with easy plane anisotropy, i.e. magnetization M in equilibrium 
state is perpendicular to the anisotropy axis N. As shown in Fig. |^ we direct z-axis along the 
magnetization vector and x-axis along the anisotropy axis. Thereby the y, z plane is the easy 
plane for magnetization vector. Assuming that resonator filed by ferrite can be considered as 
a thin slab we have the following demagnetization factors A^^: = = 0, A^^ = 1- Moreover 
following Kittel and Lax and Button we introduce an effective demagnetizing factor 



as follows: Na = 2Ka/M'^ which is directed along the x-axis. Then in this specific Cartesian 
system of coordinates the susceptibility has the following components 

_ UJoUJm 

Xxx[^) 5 

UJ;^ — UJ^ 



where 



ujQ = 'y{Ho- 4:7iM),ujm = ^tt^M, 

OOa = iHa, Ha = SrtKjM. (12) 
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Here Hq is external constant magnetic field applied along the z-axis, i.e. along the direction 
of magnetization and Ha is effective anisotropy field directed along the x-axis. 

Since for this case of the susceptibility fi^y + fiy^ = 0, expression for (|1^) and transfor- 
mation (H) are simplifying as follows 

= /ixx + — , (13) 




(14) 



If a ferrite has the anisotropy axis directed parallel to lateral sides of the rectangular 
resonator, the transformation (|Hp maps the rectangle onto rectangle. However, if the 
anisotropy axis of ferrite is directed differently, the rectangle is mapping onto parallelo- 
gram as shown in Fig. |1| (a). In particular case of squared resonator transformation ( p!^ 
maps a square onto a rhombus as shown in Fig. |l| (b) with the vertex angle of the rhombus 



a = 2arctan./^. (15) 




Because of frequency dependence of the susceptibility components expressions for eigenvalues 
of the wave equation (||) and the vertex angle (p!5[) are complicated. So we consider a 
case for which we can neglect by frequency u in denominator of the susceptibility (ITll). 
Therefore we are to imply the inequality u'^ ^ u"^. In order to fulfill this inequality we can 
apply strong magnetic field such as jIIq u; or explore ferrites with large anisotropy field 
Ha = 2Ka/M ^ LJ. If to take a size of the resonator as 10cm we obtain that, for example, the 
first thousand of eigen frequencies does not exceed the frequency about lO^^sec^^. Therefore 
we can write inequality Hq ^ lO^Oe. If to take ferrite A^io.932C'oo.o68-^e204 we have 
4itM = 3475(75, Ha ~ 125660e. Therefore we can neglect by frequency dependence in 
denominators of the susceptibility components. The anisotropy field is extremely large in 
ferrite Ba2Co2Fei2022{'2Ka/M = 2800Oe) , so upper boundary for the frequency reaches a 
value of order 10^^. 

As a result we can write the following formulas for the susceptibility components instead 
of (O) 



Ho + Ha' 

AttM 

^^^^""^ Ho-inM' 
Correspondingly substituting these expressions into Eqs (|13D and (jl^) we have 



a ^ 2 arctan 



(18) 



\(^o-4vrM)(l + ^)- 
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One can see that if the external magnetic field Hq and the anisotropy field Ha are both 
strong we obtain from ( p!7|) and (|15|) /i~l,a~7r/2. It means that the strong magnetic field 
applied along the magnetization of ferrite (z-axis) diminishes effect of anisotropy. Thereby 
we assume that only the anisotropy field is strong to obtain 



/i ~ 1, /in 



a ~ 2 arctan 







Hn - 47lM 



(19) 
(20) 



On the one hand, formula (0) shows that eigen frequencies of filed resonator coincide with 
those of empty resonator. On the other hand formula ( ^0|) gives a remarkable possibility to 
change the vertex angle of the rhombus by external magnetic field. 

While integrable, for example, rectangular billiard, and chaotic, for example, the Sinai 
billiard, represent two extreme billiards, the rhombus with at least one angle in the form 
rmT/n{m ^ 1) belongs to pseudo integrable systems |I^^18|. Integrable billiard generally 
leads to uncorrelated energy levels (Poisson statistics) and chaotical billiard corresponds to 
the Wigner-Dyson statistics 0]. Rhombus billiards are pecular as they are pseudo integrable 
systems and for this reason their statistical properties belong to another class of universality 
p!9| . Using the boundary element method pO[ Shudo and Shimizu ||16[ have found that the 



data of nearest-neighbour level-spacing distributions are described by the Brody distribution 



A= (1+ P)a, a 



expi^— as 
' 2 + p ' 

1 + /^. 



(21) 



The Brody distribution is is a semiempirical interpolation between Poisson and Wigner- 



Dyson distributions |^ . Shudo and Shimizu compared the Brody parameter P for irrational 
angles with that for the rational angles and found that the differences observed are very small. 
However, the spectral rigidity reveals the possibility of the rationality of vertex angle being 
an important signature of the level-spacing distribution. Bogomolny et al |jT^ performed 
extensive numerical calculations for a = vr/n, n = 5, 7, ...,30 and typically 20000 energy 
levels and found that the nearest-neighbour level-spacing distribution is described by the 
semi-Poisson statistics 



Ase 



-2s 



(22) 



Biswas and Jain [|l^ numerically considered the amplitude distributions P{'4') of eigenfunc- 
tions of the 7r3 rhombus. They have found that eigen functions even-even and even-odd 
relative to the x,y axes (Fig. ^ for this rhombus P(V') have the Gaussian distribution 
in contrast to the odd-odd and odd-even eigenfunctions which are identical to those of an 



equilateral triangle. These eigenfunctions are |22 



A 



COS sm 



27r(2m 



n]x 



— cos sm 



27r(2n — m)x 
3L 





2Timy 




sin 




(cos sin) 



3L 

— 27r(m + n)x 



sm 



27my 



3L 



sm 



3L 

27r(m — n)y 



(23) 
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with the eigen energies Emn = E/^{n'^ + — mn) where = 16^^7r^/(18mL^) with 
m > 2n, n = 1, 2, 3, . . . and L is the length of triangle side. This size related to The length 
of square side I related to this length by relation L = //a/2 sin(a/2) (Fig. |l|). From equations 
a = 2n/3 and ( PO] ) we obtain that this integrable case of the anisotropic square resonator 
takes place if 

Ho = 2tiM (24) 

otherwise the resonator is non integrable one. 
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FIGURES 



FIG. 1. Mapping of rectanglar billiard described by the anisotropic Helmholtz equation^ on 
to polygon (a) (transformation (^)) and of square billiard on to a rhombus (b) (transformation 

(HI))- 

FIG. 2. Schematical view of the rectangular resonator filed with ferrite where M is the mag- 
netization of ferrite and N is the anisotropy field. 
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